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Abstract 

From radio-electronics signal analysis to biological samples actuation, surface acoustic waves 
(SAW) are involved in a multitude of modern devices. Despite this versatility, SAW transducers 
developed up to date only authorize the synthesis of the most simple standing or progressive waves 
such as plane and focused waves. In particular, acoustical integrated sources able to generate 
acoustical vortices (the analogue of optical vortices) are missing. In this work, we propose a 
flexible tool based on inverse filter technique and arrays of SAW transducers enabling the synthesis 
of prescribed complex wave patterns at the surface of anisotropic media. The potential of this 
setup is illustrated by the synthesis of a 2D analog of 3D acoustical vortices, namely ’’swirling 
surface acoustic waves”. Similarly to their 3D counterpart, they appear as concentric structures of 
bright rings with a phase singularity in their center resulting in a central dark spot. Swirling SAW 
can be useful in fragile sensors whose neighborhood needs vigorous actuation, and may also serve 
as integrated transducers for acoustical vortices. Since these waves are essential to fine acoustical 
tweezing, swirling SAW may become the cornerstone of future micrometric devices for contactless 
manipulation . 
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I. INTRODUCTION 


Surface acoustic waves (SAW) have become the cornerstone of micro-electro-mechanical- 
systems (MEMS). SAW are not only useful in delay lines and convolution hlters^, but can also 
monitor temperature variations, strain^, magnetic fields^>^ and even chemical or biological 
compositional^. More recently, the growing field of microfluidics has expressed tremendous 
interest towards SAWs^*^, due to their versatility for droplet actuation^AS, atomizatioiJ^ii^, 
jettingi^ii^ or mixing!^ but also bubbles, particles and cells manipulation and sortingi^^— . 

In order to improve their actuation power, researchers have dedicated years of efforts to 
SAW focusing, with significant success^^i"— . An important side-product of this quest was the 
introduction of anisotropic zero-order Bessel functions^^. These functions are notorious in 
the fields of electronics^Z'i^S^ optics^Si^i and acoustics^Si"— for describing Bessel beams (fig. [T]), 
a large family of vortical waves spinning around a phase singularity. In this singular point, 
destructive interferences lead to the total cancellation of the beam amplitude. The resulting 
dark core is one of the most attractive feature of Bessel beams as it can trap particles^^^— . 
In the present study, we aim at generating a two-dimensional version of acoustical vortices, 
called for convenience swirling surface acoustic waves. 

In two dimensions, swirling SAW would appear as a dark spot circled by concentric bright 
rings of intense vibrations. It is tantamount to cloak the focus of surface acoustic waves, 
allowing vigorous actuation of the direct neighborhood of fragile sensors. Furthermore, SAW 
easily radiate from a piezoelectric solid to an adjacent liquid, simply by diving the transducer 
in the fluid. Swirling SAW could therefore serve as integrated acoustical tweezersi^i^S. This 
would solve one of the major shortcomings of advanced point-wise acoustical tweezers^Si^— , 
which are complex mechanical assemblies of numerous individual transducers, whereas the 
present swirling SAW generators are obtained by metal sputtering and photolithography on 
a single piezoelectric substrate. The radiation of swirling SAW in adjacent liquid might also 
be used to monitor cyclones-like flows^I in cavities by using a nonlinear effect called acoustic 
streaming. 

Two difficulties must be overcome prior to the generation of swirling SAWs. First, al¬ 
though SAW synthesis is well mastered for single transducers radiating in specific directions 
of piezoelectric materials, the design of interdigitated transducers arrays (IDTA) surround¬ 
ing a control area is much more complex. Indeed, piezoelectricity happens at the expense of 
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FIG. 1. A particular example of isotropic dark beams: the Bessel beams (eq. [T]in part |II|) with 
I = 1, kz = I and kr = 1. A) Beam cross-section with complex phase and amplitude. B) Iso¬ 
phase surfaces at W — kzZ = 0 and W — kzZ = tt in red and blue respectively. C) Iso-surface of 
3^(IFi) = —0.3 and -1-0.3 in blue and white respectively. 

isotropy, as better piezoelectric coupling requires larger substrate anisotropy^. Anisotropy 
considerably complicates the SAW propagation, leading to direction-dependent wave veloc¬ 
ity, coupling coefficient and beam stirring angle (non-collinear wave and energy vectors). 
Thanks to recent mathematical development a ^ * ^^— , SAW far field propagation is better 
handled nowadays. Nevertheless, these methods require an accurate depiction of the tar¬ 
get field in order to design the generator. The second difficulty is to define exactly what a 
swirling SAW is, especially in an anisotropic medium. Since these waves are the fragile result 
of destructive interference, extreme care must be taken in computing their propagation. 

In the present study, we use an adaptive field synthesis method in order to tackle the first 
issue. For this purpose, a sample of piezoelectric material is covered with a circular array of 
32 independent transducers actuated by a programmable electronic. Then, its vibrations are 
monitored by a Michelson interferometer. The exact input was computed by an advanced 
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calibration procedure called inverse filter— 

Getting rid of the issue of emitter design, we efficiently focused on the definition of 
swirling surface acoustic waves. Our theoretical work was essentially guided by Laude et 
ah—, who unveiled a zero-order anisotropic Bessel. In a different context (multipole ex¬ 
pansion of electro-magnetic waves for numerical computation), Filler and Martin proposed 
a comprehensive extension of Bessel functions to anisotropic medial. Our theoretical in¬ 
vestigation, described in the first part of the paper, uses the concepts of slowness surface 
and angular spectrum to fill the gap between Filler and Martin mathematical expression 
and surface acoustic waves. The next part of the paper describes our experimental setup, 
from the transducers design to the SAW measurements. The third part explains how we 
designed the emitter signal in order to synthesize swirling SAW. It provides the key steps 
of the inverse filter method adapted to the propagation of surface acoustic waves. Finally, 
a fourth section exhibits some experimental swirling surface acoustic waves. To the best 
of our knowledge, it is the first time focused surface acoustic waves are synthesized with a 
phase singularity and the associated dark spot. 


II. DEFINITION OF AN ANISOTROPIC BESSEL FUNCTION 

A classical solution to the wave equation in isotropic media is known as the Bessel beam: 

In this equation, r, 6, z, t, I, Wi, Ji, kr, kz, OJ, IF° stand respectively for the axisymmetric 
coordinates, the time, the topological order, the complex wave field value, the order Bessel 
function, the radial and axial parts of wave vector, the angular frequency and the isotropic 
swirling surface acoustic wave complex value. 

Most undergraduate courses introduce the Bessel function through series development. 
In contrast to more common waves, surface acoustic waves are the result of linear algebraic- 
partial-differential equations^, a class of problem where the use of series development, if not 
totally inefficient, would at best be cumbersome and inelegant. 

The slowness surface and angular spectrum^ constitute the basic blocks for building an 
anisotropic Bessel beam. Hence, we will first briefly review these concepts, and then use 
them to design an anisotropic swirling SAW. The main idea of these tools is to reduce the 
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problem to a superposition of plane waves. For each single direction and frequency, we solve 
the 1-D propagation eqnation, which rednces the partial differential eqnation to a set of 
ordinary differential equations, whose integration is straightforward. 

In the following, we will work at a given frequency and omit the term for clarity. 
In isotropic materials such as water, the wave speed of sound or light is independent of 
the direction of propagation. Consequently, the magnitnde of the wave-vector k = 2 tt/X 
is also a constant, and its locus versus the direction of propagation is a sphere called the 
slowness snrface. Conversely, in the case of an anisotropic material, the wave speed depends 
on the direction, and so is the wave-vector. In the reciprocal space of a 3D media, we 
call d> the azimuth and the altitude in cylindrical coordinates, so the wave-vector reads 
k{(/), Kz). The locus of this wave-vector, still called the slowness surface, then results in non- 
spherical shapes depending on the anisotropy of the material^. Bessel beams propagate 
along a specific axis z. Consequently, discussions on the surface slowness will often refer to 
kr{(t>,Hz), the projection of k on the plane normal to the propagation axis. The axial and 
radial component of the wave-vector kz and kr are linked by the direction-wise dispersion 
relation: 

kz{(pi Hz') T ^r(0) ^ (2) 

Kjz ) 

In eqnation ([2]), we wrote the coordinates in the reciprocal space Ki to distinguish them 
from ki which refer to the dispersion relation of the wave and are given physical quantities. 
For instance, can take any value whereas kz is defined only in a closed interval [kz G 
[—uj/c,-\-oj/c\ for an isotropic medinm). 

The angular spectrum is a multidimensional generalization of the Fourier transform. Since 
Fourier pioneering work, it is known that any field can be resolved into a sum of sinusoidal 
functions. The angular spectrum is a recursive application of the Fonrier transform over all 
the dimensions of the media: 


f{x,y,z) = 

r+oo r+oo r+oo 
' —oo J —oo J —oo 


F{k^, Ky, Kz)F'^^^dK^F'^«ydKyF'^^^dKz 


( 3 ) 


We can re-arrange the terms in the exponential in order to get exp{i{KxX + Kyy + KzZ)) 
snch that equation (|3D can be interpreted as a sum of plane waves. This means that any 
physical field in the media at a given frequency can be seen as a combination of plane 
waves and therefore must satisfy the dispersion relation, or equivalently lie on the slowness 
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surface. In this regard, the slowness snrface provides a frame for the wave landscape, and 
choosing the angular spectrum F{k,x, Hy, k,z) amounts to applying the color (complex phase 
and amplitude) on this frame. 

If we express the previous angular spectrum not in Cartesian coordinates but in cylindrical 
ones, we get: 


f{r,e,z) = 



In this expression, the variables 0, k-z refer to the spectral domain whereas r, 9, z belong 
to the spatial one. In order to satisfy the dispersion relation, we know that F must vanish 
anywhere except on the slowness surface, so F{Kr, (j), = h{(l), Kz)S{Kr — k-z)), with kr 

the magnitude of the wave-vector in the (x, y) plane and h an arbitrary function of (f) and 
Kz- This reduces the set of waves that can be created in the media: 


f{r,e,z) = 



( 5 ) 


At a given the integral in eqnation ([51) is the prodnct of two terms: the hrst one 
^ikr(4>,K.z)rcos{4>-e) reduced to a sum of plane waves thanks to Jacobi-Anger expansion, 

while the second one hkr provides the color of each of these planes waves. 

We construct anisotropic Bessel functions by splitting the wave angular spectrum in a 
K^-independent part and extracting its coefficients. Since 0 is the azimnth, it is a periodic 
fnnction and we can expand hkr in Fourier series: hkr = We then get: 


/(r,6»,z) = 

/ _|_QO +00 /*+7r 

l=-oo 


( 6 ) 


As mentioned earlier, the integral can be interpreted as a snm over all the Kz of some 
elementary functions. In these functions, Kz appears as a parameter instead of a variable. 

In order to highlight what in this expansion may be reminiscent of a Bessel, we need to 
write the integral expression of the Bessel function: 

r>+7r 


Jl 


1 r 

(x) = — / 

^ ^ 27r7_, 


^ilr]—ix sini 




( 7 ) 
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A trivial change of variable r] = (j) — 9 — tt/2 yields: 



We combine equation ([I]) and (jH]) to get the isotropic swirling SAW: 




h+TT 


27rz^ 


^il4>+ikrr cos{4>— 


(8) 

( 9 ) 


By analogy with the isotropic equation, we define an anisotropic swirling SAW with a 
given K,z = kz as: 



The beam in equation ffTOl) corresponds to the beam introduced by Filler and Martin^ 
for solving anisotropic multipole scattering problem, which augurs such anisotropic Bessel 
beam might be extremely widespread in Nature. 

Interestingly, any wave in an anisotropic media can be written as a combination of 
anisotropic Bessel beams W/ = 




+ 00 

ai{Kz)2m^Wf{r,e,Kz)e^'^^''dKz 


( 11 ) 


/=-oo 

In the rest of the paper, we will use the inverse filter to generate anisotropic swirling 
SAW Wf on the surface of an anisotropic piezoelectric crystal. 


III. EXPERIMENTAL SETUP 

The experimental setup is designed to be as versatile as possible, in order to allow gen¬ 
erating a wide variety of waves on an area called the acoustic scene. Starting from an X-cut 
of lithium niobate crystal, 32 unidirectional interdigitated transducers (SPUDT - IDT) were 
deposited on its periphery (see figure E]) . In order to widen the range of possible acoustic 
fields, every spot on the scene should be illuminated by all the transducers. This spatial 
coverage should be as uniform as possible on the acoustic scene. Furthermore, since any 
wave can be described as a combination of plane waves, it is essential to generate waves 
from a wide span of directions. Hence, the quality of the wave field synthesis critically de¬ 
pends on the span of plane waves provided by the source array in terms of incident angle, 
which is the angular spectrum coverage. The best way to achieve such optimal coverage is 
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therefore to gather many sources from all the directions, and dispose them radially around 
a target spot which will be the acoustic scene. This spot should also be far enough away 
such that diffraction ensures a uniform spatial coverage. These notions of optimal coverage 
are detailed further in the next section and in the appendix „ 



FIG. 2. Intedigitated transducer array used for generating the surface acoustic waves. The central 
black disc (25 mm diameter) is a gold layer acting as a mirror for interferometric measurements 
and materializes the maximum extent of the acoustic scene. Vector format image (available online) 
allows to visualize the fine structure of the electrodes. 

In order to measure the wave held on the acoustic scene, we placed the sample under 
the motorized arm of a polarized Michelson interferometer (Fig. |3]). The poor reflection 
coefficient of lithium niobate was significantly increased by the deposition of a thin layer of 
gold on the acoustic scene (~ 200 nm). 

During the design of the IDTs, special care was given to the anisotropy of the lithium 
niobate substrate. Indeed, IDTs are high quality spatio-temporal resonant elements with a 
spatial period equal to the wavelength. Any deviation from the narrow resonant bandwidth 
results in very significant loss of efficiency^. We trace in figure 0] the slowness contour of 
lithium niobate measured on the gold layer at the working frequency of 12 MHz, and compare 
it to theoretical predictions^. 

The knowledge of the dispersion relation provides the wave field radiated by a single point 


A = 632.8 nm He-Ne laser 




FIG. 3. Polarized Michelson interferometer used for scanning the displacement field associated 
with surface acoustic waves. 


source^^^: 


G(r,0) c 


( 12 ) 


With a(0) the coupling coefficient between the field to measure and the electrical potential 
(obtained when solving the SAW equations^>^) , 0(0) the beam stirring angle, /i(0) = 
cos(0 — 9)k((/))/u){(/)), h" = related to the focnsing factor. The beam stirring angle 

is solution of h'(0) = 0. 

Thanks to superposition principle, we can use the Green fnnction in eq. dED to compnte 
the acoustic field radiated by our emitter arrays. The predictions are compared to exper¬ 
iments in figure O Anisotropy strongly affects the SAW propagation, as we can observe 
beam widening (A,B), focusing (C,D) and stirring (E,F) depending on the beam direction. 

A wide-band high-power multichannel FPGA (Lecoeur Electronics) powers the 32 emit¬ 
ters with tailored numerical input. The input is specific to each desired wave-field and 
designed through the inverse filter method. 
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FIG. 4. Theoretical slowness contour (rad./mm) under a very thin gold layer— versus experimen¬ 
tally measured one. 

IV. INVERSE FILTER THEORY 

Inverse filter— is a very general technique for analyzing or synthesizing complex signals 
that propagate through arbitrary linear media. This method is especially suited for pro¬ 
totyping, because given a set of independent programmable sources, it finds the optimal 
input signal to get a target wave field. When used for this purpose, it is similar to computer 
generated holography in optics^. This method was previously shown to be among the most 
accurate ones for generating acoustical vortices in isotropic mediai2i. 

The method happens in four distinct stages (see figure Ej), (1) calibration of the trans¬ 
ducers, (2 and 3) computation of the optimal input, and (4) actuation of the sound sources 
according to the optimal input. 

In the current system, we used a set of 32 emitters, and an arbitrary number of con¬ 
trol points evenly distributed on the acoustic scene. Their density is governed by Shanon 
principle: the distance between two points should not exceed A/2. In our aquisition, we 
used a step of A/6 = 30 /rm. Moving the arm of the interferometer, we were able to reach 
individually each of these measurement points. If we call e, the temporal input of emitter 
b and Sj the temporal output of control point j located on {xj,yj}, we have for any linear 
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Theoretical 


Experimental 



FIG. 5. Influence of anisotropy on the propagation of SAW generated by single electrodes. A,C,E) 
theoretical predictions, B,D,F) experimental measurements. A,B) beam widening, C,D) beam 
focusing, E,F) beam stirring. Color (available online) represents the beam relative intensity over 
the substrate, and is not indicative of the ratio of intensity between two different transducers. 


media: 


Sj — ^ ^ hij * Ci (13) 

i 

Where * refers to the convolution product, and h^j is the time-response at control point j to 
an impulse input at emitter e*. In the spectral domain, Hij = is the Fourier transform 

of the transfer function at control point j of emitter i, and includes the propagation of the 
wave in the media. Using the matrix formalism, things get even simpler: 


S) = H\E) 


(14) 


In case where the transfer matrix is square and well conditioned, it can be inverted to deter- 
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FIG. 6. Inverse filter flowchart. Inverse filter happens in four steps: (1) recording of the spatial 
impulse response (H matrix) for all transducers. (2) Transformation of the H matrix from spatial 
to spectral domain, where the response is sharper. (3) Computation of the optimal input \E) for a 
desired output jS) by pseudo-inversion of the matrix H. (4) Generation of the signal from optimal 
input \E). 

mine the optimal input \E) from a desired output \S). However, the number of independent 
sources and control points is not necessarily the same, so H is generally not square and often 
ill-conditioned. In the appendix, we explain the reasons for this ill-conditioning from the 
perspective of the angular spectrum, and provide guidelines to minimize it. 

As soon as the value of \E) has been computed, the time-dependent input is obtained by 
inverse Fourier transform and sent to the FPGA amplifier to generate the wave field. 

V. EXPERIMENTAL RESULTS 

Bessel beams draw large interest for three main reasons: they do not diffract^, they carry 
a pseudo-orbital momentumiSSi^i^ and they exhibit a dark core (for nonzero order )^>^. In 
addition to these reasons, the zero order Bessel beam is the optimal beam focusing for a 
given aperture^i. In the following section, we start synthesizing a focused surface acoustic 
wave Wg, and then some simple first order swirling SAW >V°. We will seize this opportunity 
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to show the phase singularity and the associated dark spot. The size of the dark spot can 
easily be tuned, simply by changing the topological charge /, which is done in the third 
example with 7*^ order swirling SAWs. 

The zero order focused Wg Bessel wave phase and amplitude are traced on figure [71 It 
appears that theoretical fields and experimental ones are quite similar. In practice, we had 
to restrain the voltage amplitude of our instrument to about 10% is order to get a linear 
response of the interferometer (upper bond is about 40 nm). For high actuation power, we 
estimated the displacement amplitude based on the second bright ring. When setting the 
voltage to about 50%, we achieved displacement amplitude of nearly 180 nm. 



FIG. 7. Experimental and theoretically predicted zero order focused Wq Bessel wave phase and 
amplitude. Maximum experimental displacement is 40 nm. 

Figure in represents the first order dark beam >V° phase and amplitude. A dark core of 
zero amplitude with a diameter of 50 /rm is clearly visible at the center of the vortex, and 
matches with a phase singularity. This area is contrasted by very bright concentric rings. 
Despite some blur in the experimental measurements, a good matching between theoretical 
and experimental vortices is achieved both on shape and phase. 

Swirling SAW might be useful as integrated transducers for acoustical tweezers or micro¬ 
pumps. Tuning the topological order is essential to these applications for two reasons: it 
enlarges the first bright ring of the vortex (Olver formulai^ in eq. US]) and it increases 
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FIG. 8. Experimental and theoretically predicted first order Bessel wave phase and amplitude. 
Maximum experimental displacement is 36 nm. 

the pseudo-angular momentum of the wave^. The second effect itself generates acoustic 
streaming with an azimutal flow velocity proportional to the topological charge^Ii^. 

= /(I + 0.809 X + 0(/-T3) (15) 

In this last example, we suggest a way to increase the ring radius while maintaining zero 
azimutal streaming and keep working at the resonance frequency of the electrodes. When 
two isotropic vortices of opposite charge are combined, they result in a circular stationary 
wave pattern. In the present case, we summed two seventh order contra-rotating acoustical 
vortices Wy-|-W° 7 . The resulting field, shown in figure [9] exhibits a dark core with a diameter 
of about 500 pm circled by a crown made of fourteen extrema of amplitude. 

VI. CONCLUSION 

In this report, we proposed an anisotropic SAW version of acoustical vortices, labeled 
swirling surface acoustic waves. This implied solving two difficulties: first, the generator 
had to be designed to accommodate anisotropic propagation, and second, we needed to 
define accurately what are swirling SAW. The first problem was alleviated using a pro¬ 
grammable array of transducers controlled by a 2-dimensional version of the inverse filter. 
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FIG. 9. Experimental and theoretical predictions of the combination of two seventh order vortices 
W °7 of opposite charge. Maximum experimental displacement is 25 nm. 

while the solution of the second problem confirmed earlier theoretical predictions. We syn¬ 
thesized swirling SAW of different topological charges and large magnitude of displacement. 
This successful generation provides a pathway for integrated acoustical vortex generators 
on anisotropic substrates. Furthermore, since these beams are expected to radiate in any 
adjacent fluid, photolithography fabricated swirling SAW transducers may offer a credible 
alternative to current complicated acoustical tweezer devices made of mechanical assemblies 
of individual transducers. Beyond the specificities of acoustics, Bessel functions are very 
widespread in Nature and anisotropic Bessels may offer new analytical solutions for a broad 
class of linear anisotropic problems. 





Appendix: Optimal conditioning of inverse filtering 

Inverse filtering is a very versatile method to synthesize an optimal target field from a 
given number of transducers. As stated in section ITVl the method is not exempt of poor¬ 
conditioning which would result in large errors in the synthesized field, but some guidelines 
can significantly improve the quality of field synthesis. The poor conditioning of inverse 
filter has two roots: (1) spectral outliers and (2) redundant sources. 
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1. Spectral outliers 


At a given frequency, the wave field must fulfill the dispersion relation, which is to have 
its angular spectrum lying on its slowness surface. When the acoustic scene is a surface (in 
2D) or a volume (in 3D), this condition exactly happens. However, experimentally, there 
is always some noise introduced in the impulse response matrix, making it full rank (any 
spatial frequency can be created provided there is enough input power). Consequently, a 
first regularization is to remove the spectral outliers by sampling the target field not on a 
spatial manifold but on a spectral one, and along the slowness surface. 

Nevertheless, if the acoustic scene is a line (in 2D) or a surface (in 3D), the spectral 
condition is relaxed. Indeed, the angular spectrum of the target field is only partially known 
due to the projection of the field along the line or the surface. In 3D for instance, if the 
synthesis happens on an {x,y} plane, the system knows the values of and ky but ignores 
the ones of k^ which can then be freely chosen as long as the dispersion relation is fulfilled. 
In an isotropic medium, this results in kz = — kj^ — ky^. Note that in any case, 

> kx^ + ky^ which is the diffraction limit. Hence, spectral outliers in this synthesis 
appear beyond the A/2 boundary. 

A third example of spectral outliers is provided by piezoelectric generation on monocrys¬ 
tals. These substrates often exhibit a direction where the piezoelectric coupling coefficient 
sharply drops to zero. When this happens, no acoustic waves can be generated from this 
orientation and the associated angular spectrum coverage is barely zero. Once again, sam¬ 
pling the signal in the spectral space and excluding the zero-coupling directions avoids these 
outliers and allows an accurate synthesis. 


2. Redundant sources 

In practice, many transducers are used to ensure an efficient spectral coverage. Above 
this threshold, adding even more actuators may result in poorer synthesis quality^. Indeed, 
from the inverse filter perspective, sound sources act like a family of vectors to combine in 
order to build a target field. When two sources are redundant, the inversion operator can 
take any linear combination of them, and this indetermination will be solved by comparing 
the measurement noise associated with each source. A smart way is therefore to regularize 
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the reduced impulse response matrix obtained after removing the spectral outliers. The 
regularization can be achieved by a singular value decomposition. If two transducers are 
redundant, they will split in a singular value very close to zero and another one much more 
regular. Knowing the signal to noise ratio, it is then possible to discriminate which singular 
values originate from noise and which one are not^. 
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